Chromodynamic Weibel instabilities in relativistic nuclear collisions 
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Employing a previously derived formulation, and extending the treatment from purely transverse 
modes to wave vectors having a longitudinal component, we discuss the prospects for the occurrence 
of Weibel-type color-current filamentation in high-energy nuclear collisions. Numerical solutions of 
the dispersion equation for a number of scenarios relevant to RHIC and LHC suggest that modes 
with (predominantly transverse) wave numbers of several hundred MeV may become moderately 
agitated during the early collision stage. The emergence of filamentation helps to speed up the 
equilibration of the parton plasma and it may lead to non-statistical azimuthal patterns in the 
hadron final state. 

PACS numbers: 12.38.-t 24.85. +p 25.75.-q, 52.35.-g 
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II. FORMAL FRAMEWORK 
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In the exploration of high-energy nuclear collisions [1] , 
the possible occurrence of collective phenomena is of par- 
ticular interest. The present note reports on studies of 
color filamentation during the early collision stage. 

Within the framework of electrodynamics, Weibel [2] 
pointed out that self-excited transverse modes exist in 
plasmas with anisotropic momentum distributions and 
he derived their rate of growth from the linear response 
of the collisionless Boltzmann transport equation (also 
known as the Vlasov equation) . This treatment was later 
adapted to counterstreaming fluids of nucleons [3, 4] or 
partons [5-9] and it was found that also these systems 
possess the Weibel filamentation instability. Although 
idealized counterstreaming may not be realized, it is ex- 
pected that a significant degree of local momentum-space 
anisotropy will occur at the early stages of relativistic nu- 
clear collisions at RHIC or LHC and investigations have 
been carried out of the associated Weibel filamentation 
modes in chromodynamic plasmas [10-12]. 

The present work employs the approach formulated 
in these latter investigations in an attempt to achieve 
a more complete and quantitative understanding of the 
Weibel filamentation phenomenon in high-energy nuclear 
collisions. After recalling the most important develop- 
ments made in Refs. [10-13], we first discuss the general 
features of the filamentation phenomenon and consider- 
ing, for the first time, modes that are not simply transver- 
sally aligned. For plasmas relevant to RHIC or LHC, we 
present numerical results for the resulting growth rates. 
Furthermore, for suitable idealized dynamical scenarios, 
we extract the collective amplification coefficients and 
elucidate the importance of the rapid longitudinal ex- 
pansion as well as the equilibration caused by collisions 
among the partons. We conclude by discussing the pos- 
sible dynamical implications of these collective modes. 



The present study is based on the developments made 
previously in Refs. [10-13] and we recall here the most 
relevant elements of those. The treatment is made within 
the semi-classical transport framework in which the par- 
tons are described by their phase-space densities. An 
early review of quark-gluon transport theory was given 
by Elze and Heinz [14], while a very recent review may 
be found in Ref. [15]. 

In the present problem, we perturb a plasma of quarks, 
antiquarks, and gluons whose phase-space densities are 
uniform in space and stationary in time. Furthermore, 
they are distributed equally among the various color 
channels so the corresponding background densities are 
color singlets. For simplicity, the partons are assumed to 
be massless so their energies are E p = \p\. 

There are, in principle, separate phase-space densities 
for each quark flavor and each spin component, but since 
they all contribute additively we may consider just one 
generic quark density with a four-fold flavor-spin degen- 
eracy, Q(x,p), where x = (t,r) denotes the four-position 
and p is the quark three-momentum. Similarly, we con- 
sider just one common phase-space density for the anti- 
quarks, Q(x,p). 

In addition to being 4x4 Dirac tensors, the phase- 
space densities Q and Q are N x N color matrices in 
the SU(iV) gauge group that have singlet (colorless) and 
multiplet (colored) parts, 



Q{x,p) 
Q(x,p) 



Qo(p) 
Qo(p) 



5Q(x,p) 
5Q(x,p) 



(1) 
(2) 



The singlet parts represent the q and q background phase- 
space densities, Qo = Ifq and Qo = Ifq (where / is the 
N x N unit matrix in color space), with the phase-space 
number densities for each color being 

f q (p) = ±Tr[Q(x >P )} , f q (p) = ±Tr[Q(x,p)] . (3) 

The induced disturbances, 5Q and 5Q, which represent 
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deviations from the color neutrality, are assumed to be 
much smaller than the colorless background terms. They 
may be expanded on the N 2 -l SXJ(N) group generators 
{t a }, which satisfy Ti\t a t b ] = j5 a b (the trace is with re- 
spect to the color indices), and the individual color com- 
ponents, Q a and Q a , can be extracted by projection, 

Q a (x,p) = 2Tr[t a Q(x lP )} , Q a (x,p) = 2Tr[t a Q{x,p)} . 

(4) 

Because the group generators are traceless, only the col- 
orless part of Q, i.e. Qo, contributes to f q , while only the 
colored part 6Q contributes to Q a . 

The gluon phase-space density Q(r,p) is an (TV 2 — 1) x 
(N 2 — 1) matrix in the adjoint representation which is 
spanned by the N 2 -l matrices {T a }, whose elements are 
the SU(iV) structure constants, {T a ) bc = —if a bc- The 
identity f abc fdbc = NS ad implies that Tr[T a T b ] = NS ab . 
We may thus write 

g(x, P ) = g (p) + 5g(x, P ) = f g (p)i + T a g a (x, P ) , (5) 

where the individual elements are 

g a b(x,p) = f g {p)s ab - ifab C g c {x,p) ■ (6) 

The color singlet part, C/o, is the uniform and stationary 
gluon background distribution and f g (p) is the associ- 
ated phase-space number density in each of the -/V 2 — 1 
different color channels labeled by a. The various color 
components can be extracted by trace operations, 

f g (p) = 1 ^l^[g{x,p)] = -L-g aa ( x ,p) , (7) 

g a (x,p) = ^Tr[T a g(x,p)} = ^f abc g bc (r,p) . (8) 

The disturbances in the parton phase-space densities 
contribute additively to the induced current densities j a , 

X(r,t) = \gf ^ ^r[Q a (r,p,t)-Q a (r,p,t) 

+2i.fabcg c b(r,p,t)} , (9) 

where g is the QCD coupling constant, g 2 = AitTica 8 . 

In the collisionlcss idealization, the parton phase-space 
densities satisfy Vlasov-type transport equations, 

p»D^Q(x,p) + ^p^F^x),^^} = 0,(10) 

p^D^Q(x,p)-^p^F^(x),^^^ = 0,(11) 

p^V^{x,p) + 9 -^^ v {x),^^ = 0,(12) 

where {•,•} denotes the anticommutator. The covariant 
derivatives £> M and act as follows, 

= Id^-iglA^x), ■],V ll = Xd^-iglA^x), ■ } , (13) 



with and A^ being the potentials in the fundamental 
and adjoint representations, respectively, 

A»(x) = A%(x)T a , A»Jx) = -ifabcA^x) . (14) 

The stress tensor in the fundamental representation is 
Ffj, v = df_.Au - dyA^ - ig[Afj,, A v ], while denotes the 
field strength tensor in the adjoint representation. 

With the phase-space densities of the form (1,2,5), the 
transport equations, which are linearized with respect to 
5Q, SQ and 5g, can be explicitly solved in terms of the 
free Green's functions. Substituting these solutions into 
the expression (9) for the current, one finds the Fourier 
transform of the induced current, 



j£{k) = n^(fc) A a , v {k) , 



(15) 



where k = (w,k). The polarization tensor IP" (A;) docs 
not carry a color index, since the potential in a given 
color channel, A%, can induce color currents in only that 
channel a. It is given by 



W v {k) = g 



2 / d3 P 



(2^) 3 E p 



VUX 



p v k 



p a k a + ie 



dfcsjp) 
dp x 



(16) 

with the effective background phase-space density being 



/eff(p) = \fq{p) 



\h{p) + Nf g {p) . 



(17) 



The result (15) is similar to the electromagnetic case, 
the only difference being the replacement of the back- 
ground phase-space density of charge carriers, f e +s{p), 
by the effective background phase-space density con- 
tributing to each of the N 2 — 1 individual color channels 
a. In order to simplify the analysis, we assume that all 
the parton species (quarks, antiquarks, and gluons) have 
the same momentum profile, <ft(p), with fd 3 p(j)(p) = 1. 
So, f e s(p) ~ Pcs4>(p) where the effective parton density 

IS Pcff = \{Pq+ Pq) + Np g . 

Insertion of the expression (15) for the induced current 
into the field equation of motion, D fl F' ll/ (x) = j v (x), 
then yields an algebraic equation for A£(k), 



[k 2 g» v ~Wk v -W v {k)\A a , v {k) = , 



(18) 



which in turn leads to the dispersion equation for the 
collective modes, 



det [k 2 g^ -k"k v -W(k)] = 0. 



(19) 



However, due to the transversality of the polarization 
tensor, k f JI fJ ' v (k) = 0, the above equation involving a 
4x4 determinant can be reduced to an equation that 
involves only a 3 x 3 determinant. This simplification can 
be most easily understood in the Coulomb gauge, where 
the chromoelectric three- vector is given by E a = dtA a . 
Then Eq. (18) for A% is immediately transformed into 
an equation for E a and the dispersion equation for the 
collective modes then readily follows, 



det[k 2 (f J - k l k j 



■ u) e 



2Ji 



(*)] 







(20) 
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where the 3x3 color permittivity tensor e has the el- 
ements e % i(k) = 5 l i — w _2 Il JJ (k). Although the disper- 
sion equation (20) has been derived here in the Coulomb 
gauge, it is gauge independent since the position of a pole 
in the gluon propagator is the same for all gauges, even 
though the gluon propagator itself is a gauge dependent 
quantity (see Ref. [15], for example). 

In the subsequent presentation, we study Eq. (20) with 
the polarization tensor derived in the collisionlcss limit, 
in which the treatment is simplest. A general formal 
justification for adopting this scenario as a first approxi- 
mation lies in the fact that the frequency of the plasma 
waves is of order g, whereas that of the binary collisions 
is of order g 2 or g 4 , depending on the momentum trans- 
fer. Thus, for sufficiently small coupling constants g, the 
collective modes have a shorter characterisitc time than 
the collisions so, consequently, the collisions do not influ- 
ence the waves for short periods of time. However, in the 
present case, where the focus is on unstable modes, this 
argument may not hold, since the characteristic growth 
rates tend to be significantly slower than the frequencies 
of the stable modes. Furthermore, the coupling constant 
is actually not that small and the parton density is ini- 
tially relatively high. Therefore, the dispersion equation 
(20) derived within the collisionless idealization should be 
regarded merely as a starting point for obtaining a rough 
idea of the quantitative importance of the phenomenon. 
If it appears warranted, more refined treatments would 
be called for. 

III. GENERAL FEATURES 

In the present study, we consider momentum profiles 
4>{p) that are invariant with respect to rotations around 
the z axis as well as reflections in the xy plane (it then 
follows that 4>{~p) = 4>{p))- Then the profile function 
depends only on the magnitude of the momentum com- 
ponent along the symmetry axis, p\\ — \p z \, and the mag- 
nitude of the transverse momentum, p± = (p 2 + p y ) 1 ^ 2 - 
Due to this symmetry, we may generally assume that 
the wave vector is of the form k = (k x ,0, k z ), with 
k x = k± > and k z = fc|| > 0. Finally, since we are 
seeking modes for which the frequency is purely imagi- 
nary, we write w = i'y where 7 is real. 

The elements of the polarization tensor II(fc) are 

n« W = 9 V,(-»-| J + ^^^), (21, 

where (•) denotes the average over the profile function 
4>(p). Since H xy and Tl yz vanish in the present case, the 
determinant (20) factorizes into a y part, k 2 + 7 2 + H yy , 
which has no real roots, and an xz part presenting the 
dispersion equation of interest, £>k(7) = 0, where 

£> k = (k 2 +^+tt xx )(k 2 x+1 2 +tt zz )-(k x k z -n xz ) 2 . (22) 

Since depends on 7 only through its square, the roots 
come in opposite pairs, Wk = ±«7k- Furthermore, it can 



be seen that Z?k vanishes for 7 = since we then have 

/ n xx n xz \ = / k z (v x [---}} -k x (v x [- ■ ■]} \ 
\U ZX IL ZZ ) \k z (v z [..-]) -k x (v z [..-]) J ' 

where [•••] = [g 2 p e s{-v z d / 'dp x + v x d/dp z )/k ■ v]. The 7 
dependence of II notwithstanding, (7) appears quali- 
tatively as a fouth-order polynomial in 7, with ~ 7 4 
for large \-f\. Consequently, in addition to the double 
root at zero, there is at most one other pair of real roots, 
±7k- Thus, for any given wave vector k, the dispersion 
equation (22) has at most one positive root, 7k- 

When the wave vector is perpendicular to the axis of 
symmetry, k = (k±, 0, 0), the polarization tensor is fully 
diagonal, II IZ = 0. Furthermore, H xx > 0, so the disper- 
sion equation (22) reduces to 

= k 2 ± + 7 2 + n zz , or k\ = uj 2 e zz . (24) 

This equation has a positive root as long as k± < fc max , 
the maximum being determined by fc^ax = — ^-zz{l = 0), 

This relation shows that there are Weibel filamentation 
modes as long as the momentum profile falls off "more 
rapidly" in the transverse direction than longitudinally. 

It is readily seen that the polarization tensor remains 
unaffected under a rcscaling of its four-vector argument 
k = (w,k), Il(afc) = n(fc). Thus, since Tl(k) is pro- 
portional to g 2 p e s the entire dispersion equation remains 
unchanged if at the same time that quantity is multiplied 
by a 2 . Consequently, i/for any given value of g 2 p e ff, the 
mode with the wave vector k has the growth rate 7k, 
then in a plasma where the value of g 2 p c s is changed by 
the factor a 2 , the mode with the wave vector ak has the 
growth rate 07k- This general scaling property is very 
useful, as it makes it possible to significantly extend the 
utility of results obtained for specific scenarios. 

As k± is varied from zero to /c max (keeping fc|| = 0), 
the growth rate 7k starts out linearly from zero, has a 
broad maximum, and then decreases towards zero again. 
The initial slope, (d'j/dk±)o, is determined by 

= Vx i } ' (26) 

and it increases with the degree of anisotropy. 

For any given value of k± for which the transverse dis- 
persion equation (24) has a solution for k = (k±,0,0) 
(i.e. for k± < k max ), there exists a range of values 
fc|| < k™ ax (k±) for which the general dispersion equation 
(22) has solutions for k = (fcj.,0, fc||). [For fc|| = the 
determinant considered as a function of 7 2 , starts 
out from zero and exhibits a negative minimum before 
starting its steady growth towards ~ 7 2 . The value of 7 
for which Z?k becomes positive is the associated growth 
rate 7k- As fc|| is increased, the minimum moves inwards 
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(as does the crossing point 7k) and it eventually merges 
with the maximum at zero. We then have d 2 D^/d'j 2 = 
and k\\ has attained its maximum value fcy lax (fc_L).] 

As the direction of the wave vector k is moved out of 
the transverse (xy) plane, the direction of the electric 
field vector E starts to deviate from the z axis, as does 
the direction of the induced current j. In the non-Abelian 
plasma, the specific behavior is gauge dependent and we 
shall illustrate it in the Coulomb gauge (Sect. VI). 



IV. CALCULATIONAL SCENARIOS 





Longitudinal momentum (GeV) 



In order to have a concrete framework for the presenta- 
tion of the numerical results, we shall adopt the following 
standard values for the RHIC and LHC scenarios [16, 17], 



FIG. 1: Contour plot of the momentum profile function 
4>(p) for either the Gaussian (top) or the pQCD (bottom) 
parametrization, using a± = 300 MeV/c and an — 1 GeV/c. 
The presentation is logarithmic with four bands per decade. 



RHIC: a s = 0.3 , p cS (t = 0.4 fm/c) = 6 fm~ 3 , (27) 
LHC: a s = 0.1 , p cS (t = 0.3 fm/c) = 50 fm~ 3 . (28) 



TRANSVERSE MODES 



Obviously, they represent very rough estimates only and 
we have tried to err on the conservative side. Larger val- 
ues for either the coupling constant or the parton den- 
sity would only enhance the filamentation effect and the 
simple scaling properties of the dispersion equation (22) 
makes it relatively easy to determine what would result 
if these input values were changed. 

In our studies, we shall employ two different analytical 
parametrizations of the momentum profile function <p(p). 
The first form is a simple Gaussian, 



/Gauss 



1 



1 



27rcr z 



(29) 



al and erjf = ($,) = a 2 , 



where a\ = (p 2 x + p 2 y ) - „ x -,- „ y ^ 
Since the results will depend only on the anisotropy 
Oz/ctz, we shall adopt a fixed value for the transverse 
variance, a± = 300 MeV/c, and then vary the longitudi- 
nal width CT||. We note that (d/dp 2 ) = (d/dp 2 ) = -crj 2 
while (d/dp 2 ) = —\o^ 2 for the Gaussian profile. 

As an alternative, we employ a pQCD-motivated pro- 
file, which has a polynomial transverse falloff, 



0pQ CD ( P ) 



l 



7T (p 



2 \3 



2-KO z 



(30) 



We still take the transverse width to be a± — 300 MeV/c 
and have (d/dp 2 ) = —\v^ 2 while now (d/dp 2 ) = 

(d/dp 2 ) — — 2crJ 2 . [We have also considered replacing 
the longitudinal profile by a Gaussian in the rapidity 
Y = h \ n [(E +Pz)/(E —p z )], but such a form is unsuitable 
because of its singular behavior at p± = 0.] 

Figure 1 illustrates shape of these two test profiles. 
Generally, for similar values of a± and an , the pQCD 
form exhibits a slower transverse falloff which (as we shall 
see) enhances the degree of instability. 



In order to illustrate the typical appearance of the dis- 
persion relation, we show in Fig. 2 the growth rate 7k for 
purely transverse modes k = (fej_,0, 0) as a function k±, 
for a few selected values of the momentum anisotropy in 
matter corresponding to the initial RHIC scenario. The 
curve 7(fc_i_) rises linearly from k± = with an initial 
slope that increases with the anisotropy (see Eq. (26)). 
The curve then exhibits a maximum 70 at the optimal 
wave number k± = ko, and then decreases roughly lin- 
early towards zero as k± approaches k± = fc ma x, with 
a slope that decreases as the anisotropy increases. Al- 
though the two different momentum profiles yield quali- 
tatively similar dispersion relations, the slower transverse 
falloff of the pQCD form generally increases the degree 
of instability. As a consequence, the maximum wave 
number fc max becomes significantly larger and, perhaps 
more importantly, the maximum growth rate 70 is also 
increased, although this effect is only rather modest. 

In order to provide a more global impression of these 
features, we first consider the dependence of the maxi- 
mum wave number & max on the various input values (see 
Fig. 3). For the Gaussian profile we have 



Gauss : k„ 



9 PeE{VzPz) -? n 



1 



1 



n 7 



a- 



(31) 



which shows that there are unstable modes as long as 
&z > o x . The limiting case, fc max = 0, corresponds to 
isotropy, a 2 
file we find 



c 2 (=|(300 MeV) 2 ). For the pQCD pro- 



pQCD: fc 2 



g PcS(VzPz 



rr-. 



(32) 



Thus, evidently, /c max is now larger and the limiting case, 
fc max = 0, is oblate, ct 2 = \a 2 x (=|(300 MeV) 2 ). The 
resulting values of fc ma x are shown in Fig. 3. Since fc^^ ~ 
a s p e s, the results for the LHC and RHIC scenarios differ 



■5 



0.25 




1.0 2.0 3.0 

Wave number k T (GeV) 

FIG. 2: The dispersion relation for four initial RHIC scenarios 
(a s =0.3, p c a—6 fm~ 3 ) with various momentum anisotropies: 
the growth rate 7 k for purely transverse modes, k = (k±, 0, 0), 
as a function of the transverse wave number k± , for the Gaus- 
sian (top) or pQCD (bottom) momentum profile <f>(p). 
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Longitudinal momentum width a,, (GeV) 



FIG. 3: The maximum wave number fc max at the initial time 
io as a function of the longitudinal momentum dispersion try 
for Gaussian and pQCD momentum profiles <f>(p) in the stan- 
dard RHIC and LHC scenarios given in Eqs. (27-28). The 
approximate values given by Eq. (33) are shown for the RHIC 
scenario (dots). 



by a factor of [(ap)LHc/(c ; p)RHic] 1 ^ 2 = §, in reflection of 
the general scaling properties of the dispersio n eq uation 
(24). Furthermore, since (v z p z ) w (p z ) = y // 2/ira z , we 
obtain the following fairly accurate approximation, 



k max ~ 9 2 peS\j ~C z 
7T 



1 or 2 



erf 



1 
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(33) 



012345678 

Longitudinal momentum width a {l (GeV) 

FIG. 4: The wave number of the fastest-growing mode, ko 
(top), and the associated maximum growth rate, 70 (bot- 
tom), as functions of the longitudinal momentum disper- 
sion cr || for RHIC (solid) and LHC (dashed) scenarios, un- 
sing either Gaussian or pQCD momentum profiles having 
a ± = 0.3 GeV/c. 



More important than /c max is ko, the wave number of 
the fastest-growing modes, since these modes will become 
predominant in the course of time. It is shown in Fig. 4 
(top). After a sharp initial rise, fc exhibits a very slow 
decrease with increasing anisotropy. For a wide range of 
tTii values, the wave number of the fastest growing modes 
obtained with the Gaussian profile is typically around 
500 MeV at RHIC and 900 MeV at LHC, corresponding 
to wave lengths Ao of 2.5 fm and 1.4 fm, respectively. 
For the pQCD profile these wave numbers are typically 
100 MeV higher. We note that unless the anisotropy is 
small, we have ko fc max . 

The behavior of the maximum growth rate 70 is shown 
in Fig. 4 (bottom) and it is qualitatively similar to that 
of ko'- it exhibits a rapid initial growth followed by a 
slow decrease. The maximum growth rate is obtained 
for <7|| w 1.5 — 2.0 GeV/c and the corresponding over- 
all shortest growth times To = I/70 are 0.90 fm/c and 
0.55 fm/c for the Gaussian profile, respectively, while 
they are about 20% shorter for the pQCD profile. Thus 
the pQCD profile, with its slower transverse falloff, leads 
to somewhat larger growth rates. However, the overall 
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results are very similar for the two profile forms. In Fig. 
4 as well, the scaling feature of the dispersion equation 
(24) implies that the LHC curves are simply a factor of 
| larger than the RHIC curves. It is thus fairly easy to 
infer the behavior of the dispersion relation for any other 
particular choices of the input parameters a s and p e ff- 

Our results for the transverse modes agree qualita- 
tively with the analytical results obtained in Ref. [11]. 



VI. GENERAL MODES 

As noted above, for any given purely perpendicular 
wave vector k = (k±,0, 0) for which there exists a fila- 
mentation mode (whose induced current is then perfectly 
aligned with the symmetry axis), there exists an entire 
range of misaligned modes that are also unstable, being 
characterized by the wave vector k = (k±, 0, fcy), where 
fc|| < kf**(k ± ). 

The dispersion equation for these general filamentation 
modes is given by Eq. (22). Since D k (j) is even and 
vanishes for 7 = 0, there is at most one positive root for a 
given k. The resulting maximum values of fc|| are shown 
in Fig. 5 for selected values of the anisotropy, in both 
the RHIC and LHC scenarios. These curves delineate 
the respective spinodal boundaries in k space and the 
isotropic metric employed in the display ensures that the 
directional information is meaningful. 

For a general mode, whose wave vector k is not per- 
pendicular to the z axis, the situation is more compli- 
cated and the associated fields and currents are no longer 
simply related to k. In order to illutrate the increased 
complexity, we consider here the k dependence of the 
electric field strength E and the color current density j 
in the Coulomb gauge. (The results are gauge invariant 
for electromagnetic plasmas.) Thus, the electric field E 
is determined by the following 3x3 matrix equation, 



[k 2 (5 y - Vk? 



(k)}El(k) = 



(34) 



Since the yy term never vanishes (and the off-diagonal 
terms involving y vanish by symmetry), we must have 
E y = 0. Thus the field vector lies in the xz plane, the 
plane spanned by the wave vector and the symmetry axis 
of the momentum distribution f(p). In the aligned case, 
the wave vector is perpendicular to the symmetry axis, 
k = (k x ,0,0), and the resulting electric field is directed 
along the symmetry axis, E = (0,0, E z ). In the general 
case, the wave vector has a component parallel to the 
symmetry axis, k = (k x , 0, k z ), and thus no longer lies in 
the transverse plane. The angle formed with the symme- 
try axis, is determined by tanf?k = k x /k z , while the 
polar angle of the associated electric field vector, #e, is 
determined by 



E x 



tani?E = t=- 



k x k z - T1 X2 



kl + r 



n 2 



e z fc2 + 7 2 + n 2 



(35) 
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For small values of k z , is near 90°, while $e is near 0°. 
These directions evolve steadily as the magnitude of k z 



FIG. 5: The maximum value of the longitudinal wave number, 
fc™ ax , as a function of the magnitude of the transverse wave 
number, k± (a mode k= (k±, 0, k\\) is unstable if k\\ < fc™ ax ), 
for the initial RHIC (bottom) and LHC (top) scenarios. The 
momentum profile is Gaussian and the various values of the 
longitudinal momentum dispersion a z are indicated. The ver- 
tical and horizontal scales are equal. 



is increased, $k decreasing and -&e increasing, until they 
become equal just as the maximum value of k z is reached, 
?9 k (/c ± ,/cfp ax (fc_L)) = ^E(fc±,fc^ ax (fc_L)). It follows that 
E _L k only when k is perpendicular to the symmetry 
axis and thus it is only those modes that have a purely 
transverse character. Furthermore, there are no purely 
longitudinal modes (for which E || k). 

For large values of k x , where 7 becomes small, there 
is not much room for k z and the limiting direction does 
not deviate much from the x direction (it approaches 
the x direction when k x approaches /c max ). Obviously, 
there is most room for k z for those k x values that are 
in the region of maximum growth, while for small k x , 
where also 7 tends to zero, there is again little room for 
k z . However, since k x is now also small, the direction 
of k is significantly affected by the addition of k z . The 
resulting unstable region in k space is thus a flat disk- 
like volume (with a central depression) which is oriented 
perpendicular to the symmetry axis of the momentum 
profile (the beam axis). 

It also follows that II • k = at the spinodal boundary. 
This feature implies that the limiting direction of the 
current is perpendicular to the wave vector (and hence 
to the electric field as well). Furthermore, the direction of 
the current j turns in the same sense as the wave vector 
k, but initially at a slower rate, so their relative angle 
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never exceeds 90°, k • j > 0. 
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FIG. 6: For a lattice of wave vectors k = (k± , 0, fey ) are shown 
the direction of the wave vector, et (the upper-left line of any 
pair), together with the direction of either the electric field, 
eE (left panel), or the current density ej (right panel), in the 
standard RHIC scenario with a Gaussian momentum profile 
having cry = 1 GeV/c. These results have been obtained in 
the Coulomb gauge. 

The evolution of these various directions with the wave 
vector k is illustrated in Fig. 6. We see that in the region 
of largest amplification the degree of reorientation can 
be substantial, typically amounting to Ai? w 30°. The 
resulting modes are rather complicated to describe, since 
the directions of the k, E, and j have no simple mutual 
relationship. 



VII. AMPLIFICATION COEFFICIENTS 

The environments produced in high-energy nuclear col- 
lisions are endowed with a rapid expansion, primarily in 
the longitudinal (beam) direction z. The effective den- 
sity then decreases towards zero in the course of time, 
Pes(t) — ► 0. We shall therefore solve the dispersion 
equation at successive times t > t and thus obtain a 
time-dependent growth rate 7k (t) for each wave vector 
k. (Such an adiabatic approach is valid only for suf- 
ficiently slow evolutions and the results should there- 
fore be regarded as approximative.) Since the maximum 
wave number for which instabilities exist, fc max , is pro- 



portional to the square root of p c g, it also decreases in 
time, fc max (£)• Consequently, for each particular wave 
vector k, there is a time, tj^ lax beyond which the associ- 
ated collective frequencies iwk are real. 

Within the adiabatic approximation, the amplitude of 
a given unstable mode, C^, evolves as 



C£(t ) exp[± / lk (t')dt'] . (36) 

Jto 



The accumulated increase of the corresponding strength 
is then governed by the amplification coefficient [18, 19], 



r k = 



Jt 



7k (t) dt 



(37) 



An elementary analysis shows that if the density falls 
off as an inverse power of time, p c g ~ t^ 13 , then Tk ~ 
k 1 ^ 2 /! 3 in the limit of soft modes. Thus the amplification 
coefficient diverges unless the falloff is at least quadratic. 

In a high-energy nuclear collision, the falloff of the den- 
sity is initially approximately inversely proportional to 
time, /? rj 1, due to the rapid longitudinal expansion, 
but it then quickens (ultimately to j3 w 3) as the trans- 
verse expansion manifests itself. Thus, if carried through, 
the adiabatic treatment would yield finite amplification 
coefficients for all modes k. For simplicity, and to avoid 
a sensitive dependence on the long-time behavior for the 
soft modes, we shall employ a simple exponential falloff, 



PeffW W Peff(t ) e- 5t/ *° = p eS {to) 



[ t 



(38) 



where St = t — to is the elapsed time. This form matches 
the initial longitudinal expansion while ensuring that 
7k (t) drops off sufficiently fast at large times to avoid 
soft divergencies. It thus allows us to evaluate the am- 
plification coefficients in a manner that is insensitive to 
the long-term behavior of the collision system. 

We first calculate Tk under the assumption that the 
momentum distribution remains constant in time while 
Peff drops off in the above exponential fashion (38). This 
is expected to provide an overestimate of the amplifica- 
tion, since both expansion and equilibration act to reduce 
the anisotropy. The resulting amplification coefficients 
are shown in Fig. 7. The curves all exhibit a steep initial 
rise followed by a gentle descent for large wave numbers. 
Since the modes with small wave numbers are subject 
to amplification for a longer time, the curves have their 
maxima shifted downwards relative to the fco values for 
the respective initial scenarios. Thus, the largest ampli- 
fication coefficients are obtained for k± w 250 GeV/c in 
the RHIC scenario and for k± w 400 GeV/c in the LHC 
scenario. Their largest values are about 0.7 and 0.9, re- 
spectively, and they are reached for cry w 2 — 4 GeV/c. 

Since the pCDQ momentum profile generally leads to 
growths rates 7k that are larger than those obtained with 
the Gaussian profile (see above), the resulting amplifica- 
tion coefficients Fk are correspondingly larger, by about 
20% in the region of maximum amplification. 
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tion of motion for the momentum distribution is simple, 



0.0 0.5 1.0 1.5 

Transverse wave number k T (GeV/c) 

FIG. 7: The amplification coefficient Fk for purely trans- 
verse Weibel modes as a function of the wave number k±, 
for the idealized case when the momentum profile remains 
frozen throughout while the density decreases exponentially 
according to Eq. (38); calculated for the RHIC scenario with 
either Gaussian (top) or pQCD (bottom) profiles having the 
specified longitudinal widths a z . 



VIII. MOMENTUM RELAXATION 

In the above analysis, we have taken account of the 
decreasing density but kept the momentum profiles con- 
stant. A more realistic treatment must take account of 
the dynamical evolution of the momentum profile as well. 
Of particular importance are the rapid longitudinal ex- 
pansion and the possibility of collisions among the par- 
tons. In order to elucidate the quantitative importance 
of these agencies, we allow the momentum distribution 
to evolve as it would if subjected to a combination of 
idealized longitudinal expansion and elastic Boltzmann 
collisions, while still assuming that the overall effective 
density behaves as in (38). (Thus, although we have ne- 
glected the possible influence of the parton-parton colli- 
sions on the dispersion relation, we do now consider their 
effect on the anisotropy in the medium.) 

If we ignore the effect of the expansion and treat the 
collisions in the relaxation-time approximation, the equa- 



m 4tP,t) 



■ r [0(p,t)-0(p)] 



(39) 



where <j>(p) is the equilibrium profile (which is isotropic) 
and t c is the relaxation time. If we are only interested 
in the momentum variances, crf(t) = Jd ppf(p(p,t), then 
Eq. (39) reduces to a set of coupled equations, 



- r tf(*)-* 2 ], 



(40) 



where a 2 is the equilibrium variance. The evolution is 
then the familiar exponential relaxation, 



+ t 2 



(41) 



We now include a longitudinal scaling expansion, 
which causes the density to decrease steadily in time, 
p(t) = p(to)t /t. [The scaling scenario is boost invariant, 
so it suffices to consider what happens in a rest frame 
at the origin, where t and z are identical to the general 
variables t and 77.] Since the collision rate is inversely 
proportional to the density, t^ 1 — pcrv, it is reasonable 
to assume that it exhibits a similar decrease, t^ 1 = C/t. 
The equations of motion (40) are then modified, 



d 

dt' 
d_ 

dt 



x 



a 2 ] 



ol = 



t t Z 



(42) 
(43) 



The last term in (43) seeks to flatten the local momentum 
distribution in response to the longitudinal stretching, 
which at the same time reduces the collision rate. 

The resulting dynamics is then more complicated. For 
simplicity, we shall assume that the sum of the mo- 
mentum variances is preserved by the Boltzmann colli- 
sions, as in elastic non-relativistic collisions. The equi- 
librium variance at a given time t is then given by 
(J 2 = \(2a 2 x + crl). The equations of motion (42-43) can 
then be rewritten on matrix form, 



d_ 

dt 



1 



C -C 
-2C 3 + 2C 



(44) 



The eigenvalues of the time-independent coupling matrix, 
A T , are determined by the following secular equation, 

-3(1 + C7)A T +3C = , (45) 

and the eigenvalues of Eq. (44) are A T = — A T /3t, i.e. 



A 



.(t) = i{l + C-T[(l + C) 2 -fC^}>0. (46) 



It then follows that the associated normal variances, er 2 -, 
fall off as inverse powers of time, 



4® =e£^=e-* = . (47) 



t 
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They represent approximately the distortion a 2 z — a 2 and 
(half) the total variance 2a 2 + a 2 , respectively. 

Since the Cartesian variances can be expressed in terms 
of the normal variances, 

a\ = Af[a 2 _ - ea\] , a 2 z = Af[a 2 + + 2ea 2 _] , (48) 

with Af- 1 = l + 2e 2 and 

^ = ^{§[l + §C + C 2 ]^-§-Ic}>0, (49) 

their time evolution can be readily obtained. It is seen 
that all the variances tend to zero for any (positive) value 
of the collision constant C, thereby making it possible to 
simultaneously achieve the continual longitudinal shrink- 
age caused by the expansion, <r z — > 0, and the approach 
to isotropy caused by the collisions, <r z — > a x . 



0.7 




Transverse wave number k (GeV/c) 

FIG. 8: The effect of Bjorken expansion and Boltzmann re- 
laxation on the amplification coefficient I\, for the RHIC case 
with o\\ (to) = 1 GeV/c in various dynamical scenarios: Frozen 
momentum distributions, as in Fig. 7 (solid); flattening of the 
momentum distribution due to the stretching of the system 
(long dashes); relaxation of the momentum distribution due 
to the collisions, assuming C — 1 so t c (to) = 0.4 fm/c (short 
dashes); and both of those agencies active (dots). 

The effects of the expansion and equilibration on the 
amplification coefficients is illustrated in Fig. 8 for Gaus- 
sian momentum profiles. In view of our previous find- 
ings, we expect that the effects are similar for pQCD 
profiles. We note that the stretching produces an only 
relatively moderate decrease (about 20%), whereas the 
collisions reduce the amplification coefficients by about a 
factor of three. The combined effect is then a reduction 
of Tk by a factor of nearly four. The collision rate has 
been somewhat arbitrarily set to C = 1, corresponding 
to an initial Boltzmann relaxation time of t c (to) = to. 
Our results show that there is a considerable sensitivity 
to this quantity which might be better estimated from 
microscopic parton cascade models [20] . 

In order to further elucidate the effect of scaling ex- 
pansion and Boltzmann relaxation, we show in Fig. 9 
the effect of including both in the approximately optimal 



RHIC and LHC scenarios, for which the initial longitu- 
dinal variance is around a\\ = 2 GeV/c (see Fig. 7). As 
already noted in connection with Fig. 8, the inclusion of 
these effects lead to an overall reduction of more than 
a factor of four, relative to the simple idealized case in 
which the momentum profiles remain frozen as the den- 
sity decreases according to Eq. (38) (shown in Fig. 7). 
The results obtained with pQCD profiles are quite sim- 
ilar, apart from the values being overall slightly larger 
(about 20% for the most amplified modes). 



1.0 




Transverse wave number k (GeV/c) 



FIG. 9: The effect of scaling expansion and Boltzmann equili- 
bration on the amplification coefficient for RHIC and LHC 
scenarios. The momentum profiles are Gaussian and have an 
initial dispersion of (Tu(to) = 2 GeV/c, for which the largest 
net amplification is obtained. The solid curves are those ob- 
tained for frozen profiles (shown in Fig. 7), while the dashed 
curves are obtained in the presence of both scaling expansion 
and Boltzmann collisions (with C — 1). 

Of course, if the Boltzmann relaxation time is in- 
creased, the effect of the collisions will decrease, and vice 
versa. Since the actual collision rate is hard to assess, it 
may only be safe to conclude that the collisions are quan- 
titatively important, while the (somewhat more compli- 
cated) effect of the expansion appears to be less crucial. 



IX. CONCLUDING REMARKS 

In the present study, we have sought to make quan- 
titative estimates of the importance of Weibel instabili- 
ties in the chromodynamic plasma created early on in a 
high-energy nuclear collision. Any quantitative calcula- 
tion must rely on specific assumptions about the plasma 
environment, including its dynamical evolution which 
is presently only rather poorly known. Therefore, our 
results are correspondingly approximate. Fortunately, 
though, the simple scaling properties of the Weibel dis- 
persion relation makes it relatively easy to infer what the 
result would have been if different input values had been 
employed. Therefore, the utility of our results extends 
beyond the specific cases presented as illustrations. 
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As a concrete framework for our discussion, we have 
adopted two standard scenarios, one appropriate for 
RHIC and the other for LHC, in terms of the coupling 
constant a s , the initial effective parton density p e e{to), 
and the corresponding starting time t - Furthermore, 
we have employed axially symmetric momentum profiles 
that are either Gaussian or pQCD-motivated (leading 
to a power falloff in the transverse direction). We have 
then kept the transverse momentum dispersion constant 
at cr_L = 300 MeV/c, while exploring the dependence on 
the (local) longitudinal spread a» (for a given density and 
profile type, the results depend only on the ratio <j\\/<j±). 

Going beyond earlier treatments of these modes, we 
have permitted the wave vector k to have a component 
along the symmetry axis as well, thus extending the con- 
siderations to modes that are not purely transverse. The 
associated polarization tensor is then no longer diagonal. 
Generally the growth rates, 7k, decrease as the parallel 
component of k is increased and the resulting region of 
instability in k space is widest in the region where the 
aligned modes have their maximal growth rates. Fur- 
thermore, the electric field of a given mode, E^, as well 
as the induced current, j k , form an angle with the sym- 
metry axis, as we illustrated in the Coulomb gauge. The 
electric field then turns towards k and finally becomes 
parallel to k at the spinodal boundary, while the cur- 
rent density turns in the same sense as k. Being initially 
aligned with the symmetry axis (and thus perpendicular 
to k), it turns at first at a slightly slower rate, so its angle 
with k becomes smaller than 90°. It then starts turning 
at a faster rate and reverts to being perpendicular to k 
at the boundary, where it is thus perpendicular to E. 

The largest growth rates 7k are obtained for modes 
whose wave numbers are perpendicular to the symmetry 
axis and they have a transverse character, Elk. Over 
a wide range of <j\\ values, these wave numbers are typ- 
ically around 500 MeV at RHIC and 900 MeV at LHC, 
corresponding to wave lengths Ao of 2.5 fm and 1.4 fm, 
respectively. The corresponding growth rates are gen- 
erally larger for the pQCD profiles, but only by about 
20% for the fastest modes. However, as the density de- 
creases in the course of time, the higher wave numbers 
are progressively disfavored, as the spinodal region con- 
tracts. Therefore, the resulting amplification coefficient, 
Tk = / Jk(t)dt, which governs the accumulated degree of 
collective growth for a given k, peaks at somewhat lower 
wave numbers, namely at k s=s 250 GeV/c in the RHIC 
scenario and at k w 400 GeV/c in the LHC scenario. In 
the idealized case when the momentum profiles <p(p) are 
kept frozen in time, the largest values of Tk are reached 
for <7|| w 2 — 4 GeV/c and amount to about 0.7 and 0.9, 
respectively, for Gaussian profiles and about 20% more 
for pQCD profiles. The inclusion of a longitudinal scal- 
ing expansion reduces these numbers only moderately by 
(by about 20%), so this dynamical complication is not so 
crucial. By contrast, the inclusion of elastic Boltzmann 
collisions among the partons leads to a significant reduc- 
tion in the degree of instability, as the associated relax- 



ation drives the momentum profile towards isotropy. For 
the adopted schematic collision rate, which corresponds 
to an initial relaxation time of t c (to) = to, the reduction 
amounts to roughly a factor of three. 

Thus, overall, we find that the degree of amplification 
of the Weibcl filamentation modes is not expected to be 
spectacular for any particular wave vector k. On the 
other hand, it appears that the effect may not be neg- 
ligible either. Furthermore, it should be kept in mind 
that there are typically a large number of such unstable 
collective modes, so their combined effect on the overall 
dynamics may be significant. 

We therefore wish to conclude by speculating about 
the possible dynamical consequences of the color fila- 
mentation phenomenon. One obvious aspect concerns 
the energy dissipation. Since the agitation of these col- 
lective modes would drain energy from the background 
system, the occurrence of color filamentation presents an 
additional agency for energy dissipation. Therefore, in 
principle, to the extent that these modes are agitated, 
one may expect a correspondingly faster equilibration of 
the parton system. 

Furthermore, since the perfect azimuthal symmetry in 
an idealized head-on collision will be spontaneously bro- 
ken by the appearance of the color currents, one may 
generally expect that the emergent filamentation pattern 
will manifest itself in the angular correlations among the 
final hadrons. In particular, a non-statistical distribu- 
tion of collective energy flow will emerge along the local 
Poynting vectors associated with each amplified filamen- 
tation mode. This expectation is qualitatively different 
from that based on the parton cascade simulations [20]. 
The breaking of the azimuthal symmetry is then caused 
by jets produced in hard parton-parton interactions and, 
consequently, the effect is carried by only a few partons 
with large transverse momenta. By contrast, due to the 
collective character of the filamentation instability, the 
azimuthal symmetry breaking will be presumably involve 
a large number of partons having relatively small trans- 
verse momenta. 

It has already been speculated that color filamenta- 
tion may have observables consequences for the elliptic 
flow [21], based on the argument that the parton trajec- 
tories tend to become concentrated within the centers of 
the filaments. The conservation of phase-space volume 
then expands the momentum distribution perpendicular 
to the filaments. A corresponding quantum-mechanical 
argument can be made on the basis of the uncertainty 
relation. 

Finally, it would appear that color filamentation might 
delay hadronization. This possibility is due to the basic 
fact that no hadronization can occur in the presence of 
color currents, since the hadronic phase is constituted en- 
tirely of colorless entities and thus unable to sustain any 
colored agitations. Thus, any collective color currents in- 
duced by the filamentation phenomenon would have to 
subside before the chromodynamic plasma could trans- 
form itself into an assembly of color singlets. It would 
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thus be of interest, in a future study, to estimate how 
quickly the induced color currents dissolve again. 

In conclusion, while it appears that color filamenta- 
tion may occur in the early parton plasma, a quantita- 
tive assessment of the significance of the phenomenon 
will require detailed dynamical treatments that are not 
yet sufficiently developed. As a step in this direction, it 
might be interesting to solve the self-consistent Vlasov 
equations (10-12) in schematic collision scenarios in or- 
der to investigate how the filamcntation modes manifest 
themselves. 
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